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SUMMARY 


Let  be  given  populations  associated  with  unknown  real 

parameters  0^,  . ..,  0^  .  The  goal  is  to  find  a  population  with  a 
sufficiently  large  parameter  in  two  stages  with  screening  out  inferior 
populations  at  the  first  stage.  Both,  the  control  and  the  non-control 
situations  are  considered  simultaneously.  Let  denote  the  class  of 

permutation  invariant  randomijed  procedures  (c|>,  ty,  6),  where  at  Stage 
1,  and  ip  decide  how  many  populations  and  then  which  ones,  respectively, 
are  selected,  and  where  at  Stage  2,  after  additional  samples  have  been 
drawn  from  the  selected  populations,  6  makes  the  final  decision.  Let 
ip*  and  <5*  denote  the  natural  decisions,  i.e.  which  are  associated 
with  the  largest  sufficient  statistics.  Under  the  assumption  of  a  common 
discrete  or  continuous  type  strongly  unimodal  exponential  family  it  is 
shown  that  with  respect  to  every  reasonable  loss  function,  procedures  of 
the  type  (0,  tp*,  &*)  form  an  essentially  complete  class  within  J)  ^  . 
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1.  Introduction. 


Let  h  ,  ...»  be  k  given  populations  associated  with  unknown 
parameters  0  ,  . . . ,  8.  €  ft  ,  where  ft  c  F  is  an  unbounded  or  bounded 

1  K 

interval.  Let  the  goal  be  to  find  a  population  with  a  sufficiently  large 
parameter  after  having  screened  out  inferior  populations  at  the  first 
stage.  If  all  (all  but  one)  populations  are  screened  out  at  Stage  1,  the 
procedure  stops  and  decides  finally  in  favor  of  none  (this  one).  Other¬ 
wise,  at  Stage  2  additional  observations  are  taken  from  all  populations 
which  have  been  selected  at  Stage  1,  and  a  final  decision  is  made  among 
them.  Throughout  this  paper  our  considerations  will  be  held  general 
enough  to  cover  both,  the  control  and  the  non-control  settings.  For 
references  of  papers  dealing  with  two-stage  procedures  of  the  type  descr¬ 
ibed  above  see  Gupta  and  Miescke  (1981). 

Assume  that  samples  fX..].  ,  and  f  Y . .  J  .  ,  can  be 

drawn  from  at  Stage  1  and  Stage  2,  respectively,  i  =  1,  ...,  k  , 

which  are  mutually  independent.  Let  the  observations  from  ik  be  real¬ 
valued  and  have  a  density  c(0^)  exp(0^x)b(x),  x  €  F  ,  0^  €  ft,  w.r.t.  p, 
the  Lebesgue  measure  on  F  or  the  counting  measure  on  7.,  resp.,  1  =  1,..., 
k  .  The  function  b(x),  x  6  F  ,  and  thus  the  underlying  exponential 
family  is  assumed  to  be  common  for  all  k  populations.  Let  Ih  =  X.^  + 
...  +  X.  and  V.  =  Y.,  +  ...  +Y.  be  the  sufficient  statistics  for 
0^  with  respect  to  the  samples  of  at  the  two  stages,  and  let  their 

densities  with  respect  to  p  be  denoted  by  f„  and  g„  ,  respectively, 

i  i 

i  =  1,  ...,  k  .  Finally,  let  =  U.  +  \\  be  the  overall  sufficient 


3 


statistic  for  6^,  i  =  1,  . ..,  k  .  For  notational  convenience,  let 
£  =  (U  ,  U^),  V  =  (V^ ,  V^)  etc.  in  the  following.  For  later 

considerations,  note  that 

(l)  fQ(u)  =  cn(0)exp(0u)bn(u),  u  6  TR  ,  0  6  S  , 

g  (v)  =  c  (0)exp(0v)b  (v),  v  €  JR  ,  0  €  fi  , 

0  m  in 

where  c  (0)  =  c(0)n,  c  (0)  =  c(0)m,  0  5  ft  ,  and  where  b  and  b  denote 
n  ’  m  *  ’  n  m 

the  n-fold  and  m-fold  convolution  of  b  with  respect  to  p  . 

At  next,  let  us  give  a  precise  definition  of  a  randomized  two-stage 
procedure . 

Definition  1.  Two-stage  procedure  (<ft,  ifr,  6). 

Stage  1;  After  having  observed  U  =  u  ,  two  decisions  have  to  be  made : 

how  many  populations  should  be  selected  and  then,  which  ones,  bet  <f> 

and  ifi  be  the  corresponding  decision  functions.  Thus  let  <|>  = 

i  =  0,  1,  ...,kj  ,  where  <(k  :  -+■  [o,  l]  is  measurable,  i  =  0,  1, 

.. .,  k  ,  and  X.  .  =  1  .  Moreover,  let  f  if/  j  s  c  (l,  ...,  k}  , 

1=0  1  s,t 

|s|  =  t,  t  2  l}  ,  where  -*•  [o,  l]  is  measurable,  s  c  f  1, 

kj,  I  s  |  =  t ,  t  2  1  ,  and  L  m  ;  1  ,  t2  1  ,  and  where  I  s '  I 

s,|s|=t  s’t 

denotes  the  size  of  a  subset  s'  of  [l,  ...»  kj  . 

If  a  decision  based  on  <p(u)  is  made  w.r.t.  0  then  the  procedure 
stops,  and  no  population  is  finally  selected.  If  a  decision  based  on 
$(u)  is  made  w.r.t.  1  then  the  procedure  stops  also,  and  a  final 
decision  is  made  based  on  ^(u)  |  i  =  1*  ...»  k}  .  In  all  other 
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cases  the  procedure  proceeds  to  Stage  2. 

Stage  2:  If  at  Stage  1,  populations  with  j  €  s  =  [i^,  •••»  itJ  > 

t2  2  ,  have  been  selected,  then,  after  having  observed  V.  =  v.  , 

X1  X1 

V.  =  v.  ,  the  final  decision  will  be  made  based  on  [6.  ~(u,  v.  ,  .... 

Lt  xt  J*5  X1 

v.  )|j  €  s  j  .  Here,  6  =  [6.  |i  €  s,  sc  {l,  ....  k},  |s|  i  2}  ,  where 
k  I  s  I 

6.  :  1R  x  ]R 1  '  -*•  [  0,  l]  is  measurable ,  i  €  s  ,  s  C  (l,  k]  , 

1  ,s 

| s |  a  2  ,  and  Z  6.  =  1  ,  | s |  2  2  . 

i6s  1,3 

Let  J)  denote  the  class  of  all  such  two-stage  procedures. 

Definition  2.  A  procedure  (<f>,  \f>,  &)  €  j)  is  called  (permutation) 
invariant ,  if  the  following  three  conditions  are  fulfilled: 

J. 

<t>  is  invariant :  For  every  i  €  [1,  . . . ,  k j  ,  u  €  TR  ,  and  permutation 
a  of  (1,  ...,  k),  4>0(u)  =  4>0<ua(1),  •••*  u0(k))  and  = 

*i(uo<l>*  •**’  ua(k)}  ’ 

ifr  is  invariant :  For  every  t  €  [l,  . ..,  k],  s  C  [l,  . ..,  kj  with 
| s |  =  t,  u  €  ,  and  permutation  a  of  (1,  ...»  k),  i|)  ,  .  (u)  = 

*s,t(u(,(i)>  •••»  ua(k))  ’  where  o(s)  =  ^  sj  . 

6  is  invariant :  For  every  s  =  [i  ,  ...,  itJ  C  (l,  ...,  kj  with  t  2  2, 
•  k  x, 

i  €  s,  ii  €  1R  ,  (v.  ,  ...,  v  )  €  )R  and  every  permutation  o  of 
X1  t 

(1,  ....  k)  with  0(8)  =  s,  «0(i)f0(8)(u,  V^,  ....  v^>  =  6i>s(u0(1), 

•••’  uo(k),  va(i1)’  *••’  vaat)]  * 

Let  ^  denote  the  class  of  all  invariant  ($,  ij>,  6)  6  2  . 

Loss  Assumptions:  Let  L(£,  s,  i)  denote  the  loss  at  £  €  ft  if  subset 
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s  £  [l,  . ..,  k]  with  js|  i  1  is  selected  at  Stage  1  and  the  final 
decision  is  made  in  favor  of  i  6  s  .  Let  this  loss  function  be  invar¬ 
iant,  i.e.  L(£,  o(s),  o(i))  =  L((eo^,  ®0()c)^»s»  *■)»  i6  s££ 

k 

£l,  ...»  k},  ,  for  every  permutation  a  of  (1,  ...»  k)  . 

k  ~ 

Assume  that  at  every  fixed  6_  6  Q  with  0^  s  0^  ,  and  s  £  {3  ,  ...  k] 


with 

0  £  js|£  k-2  ,  the  following 

four  conditions  are 

satisfied: 

(LI) 

L(6,  {2},  2) 

£  L(0_,  {l>,  1) 

(L2) 

L(0_,  s  v  {1,2},  2) 

£  L(0,  i  v  {1,2},  1) 

(L3) 

L(0,  s  «  {2},  i) 

£  L(£,  s  u  {l} ,  i) , 

i  €  s  , 

(L4) 

L(0,  s  «  {2},  2) 

£  L(0,  s  «  {l},  1), 

1*1  >  1. 

Since  we  have  made  no  assumptions  with  respect  to  the  loss  LQ(60, 

k 

say,  for  selecting  no  population  at  Stage  1  at  ££  6  ft  ,  all  reasonable 

loss  functions  in  the  control  as  well  as  in  the  non-control  setting 

should  have  the  properties  assumed  above. 

Example :  The  following  loss  function  has  been  adopted  by  Gupta  and 
Hiescke  (1982)  in  the  control  case.  LQ(£)  =  0,  L(^0,  {j},  j)  =  £.(0O  -  0..), 

j  =  1,  ...,  k,  L(6£,  s,  i)  =  cjsj  +  £,(00  -  0^),  i  €  s,  s  c  [  1,  ...,  kj  with 

js|  a  2  ,  where  0Q  6  R  is  a  given  control  value,  £•  is  non-decreasing 
with  £.(0)  =  0  ,  and  c  a  0  is  the  cost  for  every  population  which  enters 
Stage  2. 

In  this  paper  our  purpose  is  to  show  that  every  (<|>,  ip,  6)  €  Ji ^ 
can  be  equalled  or  beaten,  uniformly  in  terms  of  risk,  by  (<J>,  ip,  &*), 
and,  moreover,  under  the  additional  assumption  of  a  strongly  unimodal 
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discrete  or  continuous  type  exponential  family,  by  (4>,  6*)  ,  where 

ip *  and  6*  are  the  natural  decision  functions  which  are  defined  below. 


Definition  3.  For  every  t  6  {l,  . ..,  k},  s  c  [l,  with  |s|  =  t, 

u  €  ,  i|>*  (u)  =  | B  <u)  |  -1(0)  if  s  €  (?)  B  (u)  ,  where  B  (u)  =  £  s •  | 

S  ,T  t  r~  X  "  X  ~ ~ 

s'  c  {l,  . . . ,  k} ,  |  s'  |  =  t,  max  {u^  |  i  ?  s' }  £  min  {u.  |  j  ?  s'}}  . 

Similarly,  for  every  s  =  [i^,  •••»  with  ta2,  i?s,  u?]R  , 

(v.  ,  v.  )  €  R*  and  w.  =  u.  +  v.  ,  j  =  1,  t,  6*  (u,  v.  , 

i1  it  ij  ij  i.s  -  11 

...,  v.  )  =  |C  (w.  ,  w.  )|  1(0)  if  i  6  (?)C  (w.  ,  w.  ),  where 

\  S  X1  \  9  X1  xt 

C  (w.  ,  w.  )  =  [  i |w.  =  max  (w  |r  6  s},  j  €  sj  .  Note  that  St  is 

s  ]  r 1  ■*  i,s 

only  a  function  of  (w.  ,  w.  )  . 

il  1t 


The  results  by  Eaton  (1967)  will  play  a  fundamental  role  in  our 
considerations  and  will  be  used  repeatedly.  Instead  of  mention  it  at 
every  new  occasion  let  us  point  out  now  that  in  all  relevant  situations 
Eaton's  "property  M"  is  given.  The  argument  is  always  the  same  and  can 
be  found  e.g.  at  the  end  of  Section  2  in  Eaton's  paper.  Also,  the  specific 
loss  functions  under  concern  will  always  be  invariant  in  his  sense. 


2 .  The  results . 


The  risk  of  a  procedure  (<p,  ip,  6)  6  Jl  at  £  6  is  given  by 

k 

(2)  R(£,  (♦,♦,«))  =  Lo(e)E04>o(U)  +  ^L(0,{i},i)E0  [♦1<0)*{1}fl<U>] 

k 

+  Z  Eg  [»  (U)  1_.  »  (U)  £  L(0,s,j)6.  (U.V.  )] 

t=2-  t  •»fi1....tit}8,t  j€s  j’S  "  h  xt 


I 
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Since  for  every  (<(>,  \p,  6)  €  2) R (Q_,  (<f>,  ip,  &))  =  R((0o^, 

®a(k )  ^  *  §.  €  ft  ,  for  every  permutation  a  of  (1,  ...,  k), 

it  is  convenient  to  compare  procedures  in  in  terms  of  their  Bayes 

risks  with  respect  to  permutation  symmetric  priors.  Thus  let  0  =  (  0^  , 

. ..,  )  be  from  now  on  the  random  parameter  vector  with  any  but  fixed 

permutation  symmetric  (prior)  probability  distribution  t  defined  on 
the  Borel  sets  in  ft  It  has  to  be  assumed  now  that  L(0^  s,  i)  for 
every  fixed  i  6  s,  s  S  (l,  kj  ,  is  measurable  and  integrable 

properly.  Note  that  for  most  of  our  results  we  need  to  consider  only 
priors  t  which  have  finite  supports,  where  these  two  conditions  mentioned 
above  are  met  automatically.  The  Bayes  risk  for  (<|>,  ip,  6)  £  $  under  t 
is  given  by 

(3)  r(t ,  (<j),  ip,  6))  =f  R(6,  (<J>,  ip,  6))dx(0)  =  E  R(0  ,  (<J),  <p,  6)>  . 

ftk 

Remark  1.  Note  that  for  two  procedures  in  d ^  we  have  R((),  ,  6^)) 

s  R(0,  (<P2,  ’P2,  62>),  0_  €  ftk  ,  if  and  only  if  r(T',  ($  ,  ,  6^))  s 

r(t',  (<t>2,  ip^ ,  62>)  for  every  symmetric  prior  T'  with  support  (t') 

i  k 

=  <(fW  —  Ga(k)^  I  **  permutation  of  (1,  ...»  k)},  ££  ft 

Since  we  will  compare  procedures  in  H ^  which  are  only  different  in 
the  ’P-and  6-components,  the  natural  way  to  do  this  is  to  look  at  the 
conditional  posterior  risks,  given  U  =  u  ,  for  every  fixed  u  €  1R  . 

k 

Thus  let  u  €  3R  be  fixed,  which  in  view  of  the  invariance  of  the 
problem  can  be  assumed  to  satisfy  s  u2  s  • ••  s  u^  •  Now,  for 
(<J>,  6)  €  this  conditional  risk,  given  U  =  u  ,  is  given  by 


8 


(4)  E  {R<0,  (♦,«,6))|U  =  u J  =  <t>Q(u)E  {  Lq(0)|  IJ  =  u] 


K  K  _ _ 

+  $  (u)  Z  ,(u)E  [L(0,  { i } ,  i)(U  =  u]  +  Z  $  (u)  l - 

"  t=2  t  ~  s={il5...,it} 


[  Z  E{6 .  (u,  V.  ,  V.  )L(0,  S,  j)|u  =  u,  V.  ,  ....  V  }|U  =  uj , 

i6s  3’S  X1  Xt  X1  t 


which,  after  some  standard  computations,  is  seen  to  be  equal  to 


<j>-(u)  F  L  (0)  II  fQ  (u  )dx(0)3(u) 
T0  —  J  ,  0  -  ,6  q  —  — 


■J  V 

q  i  q 


JS.  IN- 

+  <Mu)  Z  \pf (u)  f  L(0,  {i},  i)  IT  f  (u  )dx((^)0(u 
1  £=i  ^  q=2  °q  q 


K 

+  z  a  (u)  y.  z>e  .(u)  /  z  6  (u,  v  ,  ....  v  ) 

t=2  t  ~  ssfi^...,^}  S’t  “  Rt  j6s  3*s  X1  xt 


f  l(0 ,  s,  j)  n  f.  (u  )  n  gD  (v  )dx(0)  n  dy(v  )6(u)_i  , 

nk  ~  q=l  Dq  q  *s  6r  r  "  pgs  p  “ 


where  0(u)  =  J  IT  f.  (u  )dx(0)  . 
—  ,  ,  o  Q  — 

nk  *=1  q 


Lemma  1.  Let  (<p,  ip,  5)  6  1 )  be  fixed.  Then  r(x,  ($,  ip,  6))  £ 

r(x,  (<J»,  ip,  6))  and  thus  R(£,  ($,  ij>,  6))  £  R(0_,  ($,  <p,  6)),  uniformly 
in  0_  €  ,  where  ip  ie  the  same  ae  ip  except  that  ij;^  1  =  ^  , 

x  -  1,  k  • 


Proof:  The  first  assertion  follows  from  (LI)  and  Lemma  4.1  of  Eaton 
(1967),  and  the  second  one  holds  true  in  view  of  Remark  1. 
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From  now  on  we  have  to  study  the  term  following  <j>  (u)  in  (5)  in 

more  detail.  Before  we  make  use  of  the  assumption  that  we  are  dealing 


with  an  exponential  family,  it  is  crucial  for  our  further  considerations 
to  note  that  for  every  fixed  s  =  (i  ,  •••»  i^}  with  ta  2  the  following 
holds  by  Fub ini's  Theorem. 


k 

(6)  j*  Z  6.  (ij,v.  , . . .  ,v.  )  J*  L(£,s,j)  II  ffl  (u  )  H  g  (v  )dr(9)  H  dp(v  ) 
^t  j£s  ■1,s  X1  1t  ^k  q=l  q  q  r€s  r  r  p£s  p 


k  k 

=  J*  Z  6.  (u,v.  , . .  • , v ,  )  f  L(0 ,s , j  )  IT  fA  (u  )gfi  (v  )dr(0)  H  dp(v ). 

»k  ^ ll  **  nk  i-1  i  i  P 


Remark  2.  By  (6)  we  have  inflated  the  formula  by  introducing  the  v. 's 
with  i  €  s  as  dummy  variables .  This  is  necessary  to  be  able  to  proof 
our  Lemma  3  below  with  the  help  of  Eaton's  (1967)  results.  The  nextfollowing 


Lemma  2,  however,  could  also  be  proved  without  this  trick. 


For  fixed  s  =  {i^,  ...,  i^J  with  t  £  2  ,  the  term  following 

^  (u)  in  (5)  is  seen,  after  plugging  in  the  exponential  family,  to  be 

s,t 

equal  to 


k 

(7)  J  Z  6.  (u,  v.  ,  ...,  v.  )  P  L(6,  s,  j)  H  exp(0  (u  +  v  )) 

Kk  *.  q*i  q  q  q 

k  ~  -1 
c  (0  )dt(0)  II  b  (v  )dp(v  )6(u) 
n+m  q  —  ,  m  r  r  —  ’ 


k 

where  0(u)  =  f  H  exp(0  u  )c  (0  )dt(0)  and  c  (0)  =  c(0)n+m, 
-  wn-i  9  <3  n  <3  _  n+m 

n  q_i 


66  11. 
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A  change  of  variables  w^  =  u  +  v^,  q  =  1,  . k  ,  leads  us  to 


IN. 

(8)  J’  £  6  (u,  w  -u.  -u.  )  J*  L(0,s,j)  H  exp(0  w  ) 

„k  i€s  j’S  ~  *1  X1  \  Xt  Jc  q=l  qq 


c  (0  )dt(0)  H  b  (w  -u  )dy(w  )B(u) 
n+m  q  —  ,  m  r  r  r  — 

r=l 


Lemma  2.  Let  (4>,  5)  €  2^  be  fixed.  Then  r(t,  (<J>,  ijj,  6*))  £ 

r(T,  ( 4>,  ip,  6))  and  thus  R(0,  (<t> ,  ij\  6*))s  R(0^,  (0,  ip,  6))  ,  uniformly  in 

0  €  ftk  . 


Proof:  In  view  of  Lemma  1,  we  can  assume  that  =  \p  holds.  The  measure 
k 

dx(0)  =  II  c  (0  )dx(0)  is  seen  to  be  a  permutation  invariant  o-finite 
~  q=l  n+m  q  — 

•  3c 

(c  is  continuous  )  measure  on  the  Borel  sets  in  ft  .  Thus  it  follows 

•  )( 

from  (L2)  and  Lemma  4.1  of  Eaton  (1967)  that  for  every  fixed  w  €  3R  , 


j€s  3 


JV 

<5.  (u,  w.  -u.  — ,w.  -u.  )  f  L(0  ,s ,  j )  II  exp(0  w  )dx(0) 

3,S  ~  11  X1  xt  1t  k  ”  q=l  ^  *1 


is  minimized  if  6.  = 6 *  ,  j  €  s  .  This  proves  the  first  assertion, 

l.s  ],s 

and  the  second  one  holds  true  in  view  of  Remark  1. 


Remark  3.  As  we  have  pointed  out  already  at  the  end  of  Section  1,  note 
that  the  optimal  decision  function  6?  ,  j  €  s  ,  makes  use  only  of  the 


information  contained  in  (w  ,  . . . ,  w.  ) 

il  t 


Now,  for  5.  =6*  ,  j  6  s  ,  formula  (8)  reduces  to 

j  »s  j  »s 
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From  now  on,  let  only  u6  R  with  s  s  . . .  s  u^  and 
t  €  [2,  k)  he  fixed  and  let  s  £  (l,  ...»  kj  be  variable  subject 

to  | s |  =  t  .  For  notational  convenience ,  let 


(11)  Z  (w)  =  min  /  L(£,s,j)  II  exp(0  w  )dx(0),  w  €  1R*  ,  |s|  =  t  . 

‘  16s  el  ,-l  " 


Lemma  3.  Let  s  c  [1,  . . . ,  kj  with  Is  |s|  s  k-2  be  fixed  and  let 

for  p,q  6  [1,  ....  k}\  s ,  sp  =  s  «  {pj  and  sq  =  s  «  [  q}  .  Then  for 

every  w  €  ]Rk  with  w  s  w  ,  (w)  s  £  (w). 

a  —  p  q  s  —  s  — 

r  *  n  r\ 


Proof:  Let  us  look  at 

s  £  {  1,  . . , ,  k}  with 
mentioned  already,  t 
the  Borel  sets  of  fi 
the  following. 


k 


H  .(w)  =  J*  L(£,s,j)  E  exp(0  w  )dx(0^,  say, 
s’]  k  r=l  r  r 

n 

•v  k 

|s|=|s|+l,  j  (  s,  w  (  E  .  As  we  have 
is  a  permutation  invariant  o-finite  measure  on 

k 

Moreover,  for  w  €  ®  with  w  £  w  ,  we  have 

P  q 


(i)  H  .(w)  s  H  .(w),  j  €  5  . 

—  sq,q  -  sp,3  - 

This  follows  from  Lemma  4.1  of  Eaton  (1967).  This  because  for  every 
fixed  j  6  s,  we  are  concerned  with  a  fixed  size  |s|  +1  subset 
selection  problem,  where  (L3)  assures  that  the  assumptions  concering 
the  loss  function  in  Eaton  (1967)  are  satisfied. 
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(£) 


H 

s 

q 


(w)  £  H 

>q  - 


s 


p 


>p 


This  follows  also  from  Lemma  4.1  of  Eaton  (1967).  This  time  we  are 
concerned  with  a  fixed  size  1  subset  selection  problem,  where  now  (L4) 
assures  that  the  assumptions  concerning  the  loss  function  in  Eaton  (1967) 
are  satisfied. 


Pulling  together  (i)  and  (ii)  ,  the  proof  of  Lemma  3  is  completed. 


For  the  last  step  in  our  considerations,  we  assume  from  now  on  that 
the  underlying  exponential  family  is  strongly  unimodal,  i.e.  that  the 
densities  in  this  family  are  log-concave.  This  is  equivalent  to  function 
b  being  log-concave,  as  can  be  seen  immediately.  At  this  point  let  us  recall 
that  the  measure  p  is  either  the  Lebesgue  measure  on  ]R  or  the 
counting  measure  on  Z  ,  the  set  of  integers  in  ]R  .  This  is  crucial  for 
the  proof  of  our  main  result  below,  even  in  the  case  of  m  =  n  =  1  . 


Theorem.  Assume  that  the  underlying  exponential  family  is  strongly 
unimodal  and  that  the  measure  p  is  either  the  Lebesgue  measure  on  ® 
or  the  counting  measure  on  Z.  Let  (<f>,  i>,  <5)  6  be  fixed.  Then 
r(x,  (<f>,  ip*.  <$*))  £  r(t,  (<t>,  if/,  6))  and  thus  R(0,  (<J>,  ifi*,  6*))  £ 

R(£,  (<J>,  i|),  6)),  uniformly  in  6  €  ^  . 

Proof:  In  the  continuous  as  well  as  in  the  discrete  case,  strong 

unimodality  of  a  member  of  the  exponential  family  of  densities  is 

preserved  under  its  convolutions  w.r.t.  p  .  For  details  and  references, 

see  Barndorff-Nielsen  (1978),  chap.  6.  Thus  b  (x),  xgR  or  x€Z, 

m 


-Ir¬ 


respectively,  is  log-concave.  By  Lehmann  (19 59),  p,330,  this  implies 

that  b  (w-u)  has  monotone  likelihood  ratio  in  w  with  respect  to  u  , 
m 

u,w  $  1R  or  u,w  €  Z  »  respectively. 

At  the  fixed  u  €  Rk  with  £  u2  <  ...  &  and  t€  [2,  kj, 

in  view  of  (5)-(8),  (10)  and  (11),  an  optimal  (u)  equals  0  ,  unless 

S  »  vl 

s  with  Jsj  =  t  minimizes 

k 

(12)  f  Z  (w)  FI  b  (w  -u  )dy(w  )  . 

x.  s  —  ,  m  r  r  r 

]Rk  r-2 

Here  we  are  concerned  with  a  fixed  size  t  subset  selection  problem. 

8y  Lemma  3,  J£s(w)  has  the  properties  of  the  loss  function  assumed  in 
Eaton  (1967).  Thus  the  proof  is  completed  by  an  application  of  Lemma 
4.2  in  Eaton  (1967)  and  in  view  of  Remark  1, 

Corollary  1.  Under  the  assumptions  of  the  Theorem ,  the  class  of  procedures 
<4>,  ty*,  6*)  €  constitutes  an  essentially  complete  class  in  J) ^  . 

The  proof  of  Corollary  1  as  well  as  that  of  the  first  part  of  Corollary  2 
are  obvious.  The  second  part  fallows  from  Blackwell  and  Girshick  (1954), 
sec.  8.6. 

Corollary  2.  Under  the  assumptions  of  the  Theorem,  let  ($»,  \Li,  5)  €  h l 
be  fixed.  If  (4>,  \J>,  6)  is  minimax  in  then  ($»  <{>*,  6*)  has  the 

same  property  and ,  moreover ,  both  procedures  are  minimax  in  J)  . 


Remark  4.  If  one  is  looking  for  a  Bayesian  procedure  with  respect  to  a 
synmetric  prior,  then  if  there  exists  any,  there  will  also  be  one  of  the 


type  (♦,  4>*,  6*)  6  J)  ^  .  Thus  the  problem  reduces  to  optimization  of 
<(>  which,  admittedly,  will  usually  be  still  a  difficult  task.  In  the 
normal  case  (with  unknown  means  and  a  common  known  variance)  under  a 
symmetric  product  normal  prior  a  Bayes  procedure  has  been  studied  by 
Gupta  and  Miescke  (1982)  in  the  known  control  case  with  respect  to  the 
loss  function  given  in  the  Example  in  Section  1. 
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